Gapless edge states of BF field theory and translation-symmetric Z2 spin liquids 



Gil Young Cho, 1 Yuan-Ming Lu, 1 - 2 and Joel E. Moore 1 ' 2 

1 Department of Physics, University of California, Berkeley, CA 94720 
2 Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 

(Dated: March 8, 2013) 

We study possible gapless edge states of translation-symmetric Z2 spin liquids. The gapless 
edge states emerge from dangling Majorana fermions at the boundary. We construct a series of 
mean-field Hamiltonians of Z2 spin liquids on the square lattice; these models can be obtained by 
generalization of Wen's exactly solvable plaquette model. We also study the details of the edge 
theory of these Z2 spin liquids and find their effective BF theory descriptions. The effective BF 
theories are shown to describe the crystal momenta of the ground states and their degeneracies and 
to predict the edge theories of these Z2 spin liquids. As a byproduct, we obtained a way to classify 
the BF theories reflecting the lattice symmetries. We discuss in closing three-dimensional Z2 spin 
liquids with gapless surface states on the cubic lattice. 



I. INTRODUCTION 

The recent discovery of a classification^^ for various 
topological insulators and superconductors^^ relies on 
discrete symmetries of the non-interacting fermions. For 
example, band insulators respecting time-reversal sym- 
metry can be classified in two and three spatial dimen- 
sions. In both cases there are two distinct classes of band 
insulators which cannot be smoothly connected to each 
other. For the non-trivial phase of time-reversal sym- 
metric insulators in three dimensions'^, gapless states 
of Dirac fermions emerge on its surface. The gapless- 
ness is protected as long as the time-reversal symmetry 
is respected at the boundary and the bulk gap remains 
finite. 

This is very different from t he ph ysics of the fractional 
quantum Hall effect (FQHEjpSn^, which is one of the 
most well -understood and rich examples of topological 
or( j ei li3ii4] Q n i-^g ec [g e f a f rac tional qu antum Hall 

state, there is a gapless chiral Luttinger liquicP^^ which 
can never be gapped out simply because the net chirality 
of the bulk state prevents some modes from disappear- 
ing. There is no reference to symmetry needed to explain 
the stability of the edge states in FQHE. Moreover, the 
topological order of the FQHE is intimately connected 
to the gapless edge stated The connection between the 
topological bulk theory and the edge theory is eas ily un- 
derstood from effective Chern-Simon theoryS3E3 of the 
FQHE. In this effective theory description, the edge de- 
gree of freedom is encoded in the gauge invariance of the 
Chern-Simons theory with an open boundary. 

However, Z2 spin liquids, which are also one of the 
well-established examples of topological order, do not 
have gapless edge states in general. The effective BF 
theor y 18 " 20 ! description of Z2 spin liquid does have an 
edge degree of freedom as in the Chern-Simons theory: 
the edge theory is a pair of chiral fermions propagating 
in opposite directions. In the absence of any symmetry, 
the se tw o fermion modes can backscatter to open up a 
gap 1 18 * 19 ! in contrast to the "chiral" Chern-Simons theory 
of FQHEs. This is true as far as no symmetry is imposed 



and can be generalized to any 'doubled' theory 21 . How- 
ever, it is known that the ed ge the ory of the Abelian 'dou- 
bled' theory can be gapless^ES when [/(l)-charge con- 
servation and time-reversal symmetry are present. The 
physical example of this is th e fractional (and integer) 
quantum spin Hall effectP^iHIl! As an Abelian 'dou- 
bled' theory can be re-written formally as the BF theory 
(at the level of the Lagrangian) and a Z2 spin liquid is 
described by BF theory one might think that the edge 
theory of Z2 spin liquids can be also gapless if the time- 
reversal symmetry is imposed. This is, in fact, incorrect 
for Z2 spin liquids. 

The reason why Z2 spin liquids fail to have a gapless 
edge state is traced back to the differences in the charge 
lattice of the compact gauge theory (or equivalently, the 
differences in the allowed operators of the edge theory 
due to quantization). We will see that the Z2 conserva- 
tion of charge and vortices in Z2 spin liquids is crucial. 
We will expose this 'structural' difference of the formally 
identical theories in the subsequent discussion in this pa- 
per. Nonetheless, the understanding of the gapless edge 
states of the fractional spin Hall effect indicates that the 
doubled theories can have gapless edge states. Thus, it 
implies that Z2 spin liquids will have gapless edge states if 
the correct symmetries are imposed on top of the topolog- 
ical order. In this paper, we will show that translational 
symmetry can stabilize the gaplessness of Z2 spin liquids 
in certain cases. We make a direct connection between 
the microscopic structure of the physics and the effective 
BF theory to confirm the gaplessness of the edge states. 

Another motivation of this pap er is to study the ef- 
fect of the translational symmetry^l^QI on the so-called 
"intrinsic topological order" . The intrinsic topologi- 
cal orders in gapp ed phases are featured by long-range 
entanglement!^^, fractional excitations^!, and topolog- 
ical degeneracies^. Examples of intrinsic topological or- 
ders are the fractional quantum Hall effect states^ and 
gapped quantum spin liquids^. Intrinsic topological or- 
ders are stable against any weak perturbations in the 
presence of no symmetry, to compare with "symmetry- 
protected topological orde r". T he symmetry-protected 
topological (SPT) phases^i2135| i ncm de topological in- 
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sulators and topological superconductors^^ as its out- 
standing examples. One SPT phase is a gapped phase 
with no fractional excitations and no topological degen- 
eracy, which has gapless boundary states protected by 
symmetry. Most importantly, a gapped SPT phase can 
be continuously connecte d to the trivial phase if the sym- 
metry is brokerP- 3 °l 36 l 37 [ A natural question is: what 
kind of role does symmetry play for intrinsic topolog- 
ical orders? It turns out that different symmetry en- 
riched topological (SRT) phases^^ can emerge from the 
same intrinsic topological order, such as different classes 
of fractional topological insulators 2 -^ with U(l) xi Z^ sym- 
metry. Therefore fractional topological insulator a 22 ! 38 ! 39 ! 
and other SRT phases exhibit interesting interplay of 
fractionalization and symmetry. Here, we show that by 
imposing translational symmetry new physics emerges in 
Z2 spin liquids such as gapless edge states and the Ma- 
jorana zero modes at the lattice dislocations. 

The rest of the paper is organized as follows. In section 
|nj we will study the general edge theory of Z2 spin liq- 
uids and show that the edge theory is generally gapped 
when there is no symmetry imposed on top of the Z^ 
topological order. We also compare the edge theory of 
Z2 spin liquids to the edge theory of the quantum spin 
Hall insulator. In section III, we impose translational 
symmetry on Z2 spin liquids. We reveal the structure of 
Wen's plaquette model, which has a gapless edge state, 
and generalize it to construct mean-field Hamiltonian of 
Z2 spin liquids with gapless edge states. We also study 
the edge theory in detail, including the condition for the 
existence of the gapless edge states and its stability. In 
section IV, we encode the lattice translational symme- 
try into the BF theory and make a connection between 
BF theories and the Z2 spin liquids found in section III. 
We present a generalization to a three-dimensional model 
in section V, then summarize our results and raise open 
questions in section VI. 

II. GENERAL DISCUSSION OF EDGE THEORY 
OF Z2 SPIN LIQUIDS 

In Z2 spin liquids, there are two low-energy excitations, 
namely the spinon and the vison. Spinons (visons) carry 
electric charge e ~ 1 (magnetic charge m = 1) in the 
underlying Z2 gauge theory. They can be described by 
BF theory with the charge lattice e € Z2 and m € Z2. 
We denote the electric current J M and magnetic current 

which couple to a M and minimally Then, the low- 
energy theory for the Z2 spin liquidlM^o] j g 

L = 1 e^ x a^d u b x - - b»f + 0(da, dbf . (1) 

As noted before, a formally similar BF theory emerges 
as the effective theory for the spin Hall effect with a 
Kramer pair of counter-propagating gapless edge modes. 
To manifest this similarity, we diagonalize the BF La- 
grangian into two copies of Chern-Simon theory, i.e., we 



write Afj_ = + b^ and = — b^ to obtain 
L = ~s^ x A ti d v A x -~e^B tl d v B x -A. X -B-T, (2) 

47T 47T 

where the dot product [A ■ J = A^J^) are understood. 
We introduced the source currents \ = an d t = 
for A^ and B^ . The charge lattice associated with x an d 
r can be deduced from that of J and j (i.e., of e and m). 
We begin with (e, m = 0) where e is defined modular 
2. For e = and m = 0, (x, t) = (0, 0). For e = 1 and 
rn = 0, (x,t) = (1/2, 1/2). For e = 2 - mod 2 and m = 
0, (X T ) = (1)1) ~ (0,0). Similar consideration shows 
that (x, t) = (1, —1) ~ (0,0). Thus, in the diagonalized 
doubled Chern-Simon theory, (x,r) = (1,±1) ~ (0,0) is 
the same as (e, to) = (2, ±2) ~ (0, 0) in BF theory. Now, 
we construct the edge theory from Eq.([2]). 

Ledge = ip R (dt - vd x )if) R + tp\{d t + vd x )ip L (3) 

This edge theory, in general, is unstable to opening up 
a gap 1 18 * 19 ! . This can be easily seen from Z2 conserva- 
tion of spinons and visons. The Z2 gauge charge conser- 
vation allows us to add the mass term ~ ipRipLiVLipR 
and ipR'&Li'&li&L- This is because i^r^l > ^l^r, or 
iPr4>l, V'flV'I carry charges (%, r) = (1,±1) which are 
equivalent to no charge (x, r) = (0,0) in Z2 theory. 
Hence, we conclude that there is no protected edge state 
for a general Z2 spin liquid without additional symme- 
tries, independent of microscopic constructions. How- 
ever, this could be changed when the Z2 spin liquid is 
supplemented by a symmetry. Additional symmetries on 
the topological order can restrict the form of the mass 
terms and stabilize the gapless edge state. 

We now show that the time-reversal symmetry cannot 
stabilize the 'gaplessness' of the edge states of Z2 spin 
liquids (if it does not possess any 'strong index^ for the 
underlying fermionic spinons). We take the edge theory 
Eq.(|3]) and consider the time-reversal symmetry opera- 
tion on the A^ and B^. Due to the fact that under time 
reveral T : A^ — > (-1)^ and B^ -> (-1)M M (with 
the definition (-1)^ = 1 for \i = and (-l)^ 1 = -1 for 
/i = 1,2) under the time-reversal operation T, the time- 
reversal operation effectively acts as the exchange of two 
fermionic fields, i.e., T : ijjR — > i/jl and tp^ — > e l9 ipR up to 
the U(l) phase factor e t6 '. We can take 9 = it or to be 
consistent with the time-reversal operation T 2 — (— 1) f 
or T 2 — 1, where Nf represents the total fermion number 
operator. This operation should be supplemented with 
v — > —v. Hence, the kinetic term for the edge theory 
Eq.(|3]) is time-reversal symmetric. We discuss the two 
different cases T 2 = 1 and T 2 = (—l) N f independently 
First, for T 2 = 1 the time-reversal symmetry allows mass 
terms of the form -0^-0^, with an equal amplitude for 
ip L ipR, which are capable of gapping the edge. Other 
terms such as V'flV'L are n °t allowed. As a whole, the 
time reversal symmetric edge theory for Z2 spin liquid is 

L = L e dge{lpR, Iph) + m(lp R lp L + V4,V>fi) + h.C, (4) 
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where L edge is the kinetic term Eq.Q. For T 2 = (-l) Nf , 
ip R ipL and ifj^ L ifjR are not allowed. Instead, there are al- 
lowed terms such as ipR^L, and they are enough to gap 
out the edge state. Here we assumed that there are no 
strong indices for the underlying spinons. The strong in- 
dex includes the Z2-indices for Dili class and Z- indices 
for C/D classed^ for band structures of the Schwinger 
fermion d " ' (fermionic spinons). When the strong in- 
dex of the time-reversal symmetry is nontrivial, then this 
phase automatically has the gapless edge state which is 
either a pair of helical Majorana modes (Dili class) or 
chiral Majorana modes (C or D classes). However, these 
phases are not described by BF theory-^ and will not be 
discussed in this work. 

For the quantum spin Hall effect on the other hand, 
we have £7(1) charge conservation and the time-reversal 
symmetry T with T 2 = (—l) Nf . The two conditions 
exclude^ the two mass terms ip R il>L, ipRipL which were 
allowed in Z2 theory. Hence, the edge excitations of the 
quantum spin Hall state are gapless, in contrast to that 
of Z2 spin liquids. 

In the remaining sections of this paper, we will consider 
the effects of the translational symmetries on the states 
localized at the edges of Z2 spin liquids. We find that 
translational symmetry can stabilize the gapless modes 
at the edges of some Z2 spin liquids. 



III. WEN'S PLAQUETTE MODEL AND MEAN 
FIELD THEORY OF Z 2 SPIN LIQUIDS 

Wen's plaquette modeP^ is an exactly solvable model 
on the square lattice where the spin degrees of freedom 
are fractionahzed into Majorana fermions. The model 
can be formulated as the pure Z 2 gauge theory^ 



duces to 



u, 



13 > 



(5) 



p <ij>eP 



where P denotes a plaquette of the direct lattice. Here, 



±1 is the Z2 gauge field on the link < ij >, and 



thus n<ij>ep Uij = Fp is the field strength on the pla- 
quette P. A faithful representation of this gauge theory 
is Ui^+ X = i^i,xXi+£,x, Ui : i+y = i^i,yXi+y,y with trans- 
lational symmetry. There are four Majorana fermions 
(Aj >x , \i,y, Xi,xi Xi,y) P er site, and the dimension of the 
Hilbert space per site is 2 4 / 2 = 4 states per site (see 
Figjl]). To faithfully represent a spin- 1/2 system, we need 
to halve this Hilbert space, i.e., the Hilbert space per site 
should have dimension 2. This is done by taking the Z2 
redundancy into account: the theory Eq.([5]) is invariant 
under TJ%j — > SiUijSj, Sj = ±1. The model can be more 
clearly represented if we introduce two complex fermions 



fi.u 



*Xi 



fi.d ~ ^i,V 



l Xi,yi 



(6) 



such that the spin S% on the site i is represented via 
lEa./s/i.a^A/S- Then > the Hamiltonian Eq.Q re- 



H Q ^ ^ 3i,x$i+x,y$i+x+y,x$i-\-y.y- 



(7) 



From the model Eq.([5]), we can immediately see the ex- 
istence of the edge degree of freedom. As the Majorana 
fermions in bulk will be paired within the plaquette, there 
will be dangling Majorana fermiond^ at the boundary 
because of the edge c uts a plaquette in half (see Fig[l|). 
It was already noticecP^E3 that this model (with g > 0) 
has a flat band of edge states on the boundary along 
x- and y- directional cut (see Fig[Tj). Hence, to look for 
gapless edge states for Z2 spin liquids protected by trans- 
lational symmetry, we can consider similar models which 
supports dangling Majorana fermions on the edge. We 
will generalize the structure of Wen's model to gener- 
ate a series of mean-field Hamiltonian for Z2 spin liquids 
with gapless boundary states in two and three spatial 
dimensions. We also discuss the stability of the gapless 
boundary states. 

We now look into Wen's model carefully. To illuminate 
the underlying structure of Wen's model, we study the 
mean- field Hamiltonian for the fermionic spinons (JfjJ. 

H?nean(\_fi,ui fi,d\) ^ ^ fi,T ~t~ fi.aVij fi,T 

(8) 

We notice, by plugging Eq.Q, that f^ u and f^d are de- 
coupled completely and form the one-dimensional Majo- 
rana fermion chains^ (in the weak-pairing phase) along 
x- and y- axis. Explicitly, we now have 



1 {{fi,u, fi,d}) — hiD,x({fi,u}) + hiD,y({fi,d})i (9) 



where /ii£>,ft({/}) is the one-dimensional Kitaev model 
along h in the weak pairing regime. This explains why we 
have dangling Majorana fermions on the boundary along 
x- and y- directional cut of the lattice (See Fig[I]). Wc 
will see that there is a gapless Majorana mode formed by 
the dangling Majorana fermions as long as there are odd 
number of Majorana fermions per unit cell on the edge. 
Before discussing the nature of the gapless edge states, 
we notice that we can easily generalize Wen's model by 
deforming one of /iid, «({/}) Eq.([9]) into the trivial phase 
while keeping the constraint 'one spinon per site'. Then, 
this Z2 spin liquid will break the C4 rotational sym- 
metry of square lattice while keeping the translational 
symmetry. For example, let us consider the case where 
hiD,i({fi,u}) is replaced by the trivial <2({/i,«}) (this 
Hamiltonian Q({f}) is not necessarily one-dimensional 
but is fully gapped): 



Q({h,u}) + hx D Mfi, d }) (10) 



has dangling Majorana fermions if the edge is not parallel 
to y. Similarly, we can generate a Z2 spin liquid which 
has dangling Majora na ferm ions if the edge is not parallel 
to x. For Z2A phasePEl whose spinon band structure 
preserves translation symmetry explicitly, we can have 
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(A) 



(C) 



(B) 




FIG. 1. Illustration of Wen's plaquette model and its edge 
state. (A) Plaquette operator ]~I<ij>gp^i- A plaquette op- 
erator pairs up four Majorana fermions around the square. 
Here, blue circles represent four Majorana fermions per site, 
and red lines connecting two blue circles represent Uij appear- 
ing in the plaquette operator. (B) In terms of the complex 
fermions / = A ± i\ in Eq.(6j, Wen's plaquette model can 
be understood as the stacked one-dimensional Kitaev chains 
parallel to x and y. (C) Edge states of Wen's model where the 
edge is along y direction. The plaquette operators in Hamil- 
tonian Eq.Q cover all the Majorana fermions except a single 
Majorana fermion per site on the edge. The blue circles rep- 
resent dangling Majorana fermions while red circles are Ma- 
jorana fermions covered by the plaquette operators. These 
blue dangling modes will hop around by the magnitude ~ t 
when we are away from the exactly solvable point. (D) The 
dangling Majorana fermions can be equally well understood 
via the stacked Kitaev chain picture. When the Kitaev chain 
is cut, the edge always hosts a single Majorana fermion per 
chain. It is clear that the edge states of (C) and (D) are 
exactly the same. 



another possible Z2 spin liquid where we align the Kitaev 
chain along n — x ± y (we will see later that this state 
has the same characteristics to the Wen's model) . Hence 
we have three translational symmetric Z2A spin liquids 
with the gapless edge states. Though these Z2A spin 
liquids have the same Z2 topological order and the same 
translational symmetry, we can disti nguish them further 
by looking into the crystal moment a 1 19 * 29 ^ of the ground 
states and degeneracies, and the edge states. In the next 
section, we use these crystal momenta and degeneracies 
to find the effective BF theory for the Z2 spin liquids. 

Now that we have some conditions for the existence of 
the edge states, we study their detailed features. We first 
discuss the relation between the edge degrees of freedom 
and the bulk gap. We argue that there is a gapless edge 
state as long as the bulk is gapped and the translational 
invariance at the edge remains unbroken (See Fig(2|. To 
remove the boundary Majorana fermions, we need to pair 
them up into complex fermions. Due to the translational 
symmetry, a Majorana fermion on the boundary cannot 



pair up on the boundary if there is one single Majorana 
fermion per unit cell. What the Majorana fermion can 
do, instead of dimcrizing on the boundary, is to tunnel 
through the bulk and pair-annhilate with the Majorana 
fermion on the opposite boundary. However, the tun- 
neling lengtipS is ~ A| where A is the bulk gap. 
Hence, we need to close the gap and go through the bulk 
phase transition to lose the dangling Majorana fermions 
at the boundary. In the next section, we will construct 
the effective BF theory for these phases, and the effec- 
tive theory should describe a finite range of the phase 
diagram as it is gapped and the edge state of the topo- 
logical field theory is stable as long as the gap remains 
finite. 

The above 'gaplessness' of the edge states depends on 
the number of Majorana fermions per unit cell on the 
boundary. If there are even numbers of the dangling Ma- 
jorana fermions per unit cell, then Majorana fermions 
can pair up on the boundary without breaking the trans- 
lational symmetry along the boundary. The 'parity' of 
the number of the dangling Majorana fermions on the 
arbitrary boundary can be easily computed, and there is 
an odd number^S of the dangling Majorana fermions per 
unit cell on the boundary if 



Jft = — • G = 7r mod 2-7T, 



(11) 



where n is the vector defined for the Hamiltonian of the 
Majorana fermion chain /iiD,n ({/})) and G is the recip- 
rocal vector orthogonal to the boundary. If there is more 
than one chain, we simply add up all Jn mod 2-7T. Then we 
have gapless edge modes protected by the translational 
symmetry if the sum is 7r mod 27r. In fact, the same index 
is used to study the localized states at the di slocat ions 
of weak topological insulators/superconductors^^. In- 
tuitively, h ■ G in Eq.(ll| can be thought as the "flux 
density" of Kitaev chains (which are parelle to n) passing 
through the area perpendicular to G (which is analogous 
to the definition of the flux through a surface in the el- 
ementary calculus, i.e., F — E ■ S defines the flux of the 
vector field E passing through the area perpendicular to 
the vector S). Due to the 'parity' (or even/odd- ness) of 
the dangling Majorana fermions, we need to take mod by 
27T with the appropriate normalization. To demonstrate 
this, we work out a few examples applying Eq.(ll). 



- example 1. edge state of Wen's plaquette model with 
the edge parallel to y. As the edge is along y-direction, 
the reciprocal vector defining the edge is G = (2tt, 0). In 
Wen's plaquette model, we have the two Kitaev chains 
where each chain is identified with hi — (1,0) and 
fi2 = (0, 1). We immediately have Jf ll — tt and J,~ l2 — 
and thus we have a dangling Majorana fermion per site 
because Jn t + J„ 2 mod 27r = n. This agrees with the 
previous intuitive understanding depicted in FigjT] 

- example 2. edge state of Wen's plaquette model with 
the edge parallel to x + y. This edge can be identified 
with the reciprocal vector G = (2tt, —2tt). As before, we 
have the two Kitaev chains where each chain is identified 
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with hi = (1,0) and hi = (0,1), and thus Jf ll = tt, 
Jf L2 = — tt and Jf L1 + Jn 2 mod 2ir — 0, i.e., there are 
even number of Majorana fermions per unit cell at the 
boundary and hence there is no gapless helical Majorana 
mode protected by the translational symmetry. Indeed, 
if we look at the edge defined by G — (2ir, ~2tt), there are 
two Majorana fermions (per unit cell on the boundary) 
coming from the Kitaev chains along x- direction and 
Kitaev chains along y- direction on the boundary which 
immediately implies that there is no stable gapless edge 
state. 

If we have a single Majorana fermion per unit cell, we 
can study the spectrum of the boundary modes. When 
the perturbation away from this ideal Hamiltonian is 
given, the dangling Majorana fermions will hybridize 
with the nearest neighbors and start to disperse 

Hedge = H ^ VjVj+l, 

j 

where the position of the boundary Majorana fermion is 
labeled by the index jeZ and t is the effective hopping 
parameters on the boundary. The spectrum of this Ma- 
jorana fermion is given by E(k) = 2isin(fc) for k G (0, tt) 
(this Majorana fermion problem is, in fact, related to the 
fermion doubling problem) . To gap out the spectrum, we 
need a perturbation with a matrix element connecting 
k = and k = tt. However, this interaction doubles the 
unit cell on the boundary (similar to the perturbations 
required to eliminate the surface of the weak topological 
insulato j 52 * 53 !) , and hence it is prohibited by the lattice 
translational symmetry. We note that the right-mover 
and the left-mover of the gapless edge theory have dif- 
ferent center of mass momentum. This feature is also 
reproduced by the effective BF theory. 

From this edge spectrum, the counting of Majorana 
fermions, and the translational symmetry, we now show 
that this edge theory mimics the stability of the edge 
state of a quantum spin Hall insulator against the 
coupling to ordinary one-dimensional Luttinger liquids. 
When the edge state of the spin Hall effect interacts with 
ordinary Luttinger liquids, the edge states are recon- 
structing themselves to the edge theory of spin Hall in- 
sulators with the renormalized coefficients. For the edge 
of the Z2 spin liquids, we can couple our edge theory 
to the quantum Ising chain at the critical point, which 
is a gapless helical Majorana state (see Fig(2|. The he- 
lical modes (the left-mover and the right-mover) of the 
quantum Ising chain are located at k — tt. When the 
helical modes are placed on top of the boundary Majo- 
rana modes of Z2 spin liquids, the right-mover at k = 
is intact as it cannot interact with the modes at k = tt. 
At k — tt, there are three modes (two from the criti- 
cal Ising chain and one from the original QSH edge): one 
right-mover, two left-movers. The translational-invariant 
potential will have matrix elements connecting k = tt 
to itself, allowing one right-mover and one left-mover to 
pair up. Hence, we are left with one left-mover at k = tt 
which is reconstructed from the helical Majorana modes 



of the quantum Ising chain and the original left-mover 
at k — tt. Thus, the boundary remains gapless even if it 
interacts with the critical Ising chain. It is not difficult 
to see that other systems, such as helical Dirac edge, can- 
not gap out the original edge modes as they have even 
numbers of the helical Majorana fermions. The easiest 
way to understand this behavior is as follows; the size of 
the Hilbert space per site in the quantum Ising chain is 
2, and we need two Majorana fermions per site to repre- 
sent the Hilbert space. This implies that there are three 
Majorana fermions per unit cell (see Figj2]), and thus we 
expect the spectrum to be gapless after the dimeriza- 
tion of two of them. Hence we can conclude, based on 
the counting argument and the translational symmetry, 
that the edge is robust against to coupling to any one- 
dimensional gapless system. 

It is also interesting to note that this edge spectrum 
implies that we will have a zero energy Majorana state 
at the lattice dislocation in Z2 spin liquids^. For Wen's 
plaquette model, the dislocation of the Burger's vector x, 
y will trap the zero energy Majorana states. On the other 
hand, we will have the zero energy state of the Burger's 
vector to orthogonal to h if Z2 spin liquids contains a 
Majorana fermion chain along h defined above in this 
section. 

The discussion in this section is based on the mean- 
field Hamiltonian for the Z2 spin liquids. However, the 
spectrum is fully gapped, and hence the fluctuation over 
mean-field solutions are supposed to be small. Still, one 
might think that it is unclear if our result remains the 
same when both four-fermion interactions and gauge fluc- 
tuations are included. These questions, perhaps, cannot 
be answered in the mean-field theory. Therefore we take 
a different path to answer the question. In the next sec- 
tion, we will find the effective BF theory for these spin 
liquids, and the topological BF theory should describe a 
phase instead of a point in the phase diagram. Hence the 
BF theory predicts the phase to have gapless edge states 
on the boundary protected by translational symmetry as 
long as the gap remains finite. 



IV. EFFECTIVE BF THEORY AND ITS 
CHARACTERIZATION 

In this section, we will construct the effective BF the- 
ory of the Z2 spin liqui ds con sidered in the previous sec- 
tion. In previous studies^^ZI it was shown that the trans- 
lational symmetry restricts the form of the mass into the 
specific form so that the edge state becomes gapless. And 
the same effective theory captures the crystal momenta 
of the ground states and topological degeneracies in a 
square lattice with periodic boundary condition. Here 
we extend previous results in the way that we clarify 
the spectrum of the edge theory and its connection to 
the microscopic model discussed in the previous section 
(in fact, the spectrum of the edge states is crucial for 
the gaplessness as we saw in the microscopic discussion) . 
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Furthermore, we demonstrate the way of classifying BF 
theory in the presence of lattice symmetries. 

The essence of the effective theory is to encode the non- 
trivial transformations for the gauge fields in BF theory 
(TlJ) under translations (in the Coulomb gauge ao = 60 = 

of 

tx,y ■ CLi — !• bi, hi di, (13) 

leaving the theory Eq.Q invariant. t x ,y denote lattice 
translations along x and y directions. This directly im- 
plies that the edge theory^ should be symmetric under 
ip' L — s- ipL and Vl — > "0l m the edge theory Eq.([3|. This 
symmetry, hence, restricts the form of the masses for the 
edge theory into 



(14) 



To see the effect of this particular mass term Eq.(|14|, 
we introduce tpR = r)R + i\R and ipr, = xl + if}L where 
Xr/LtVr/l are real Majorana fermions. Then, the mass 
Eq.(14| gaps out \l and one of two right-moving Ma- 
jorana modes XRiVR ( e -g- XR if a is a real number), 
and hence we are left with one right-mover and one left- 
mover. It is not difficult to see rji, is located at k = tt as 
ipL transforms to tj)' L under unit lattice translation. 

We will follow this lesson to construct and (partially) 
classify the effective BF theory for the Z2 spin liquids 
discussed in the previous section. However, the classifi- 
cation in this section is based on the educated guesses 
constructed in the previous studies. The philosophy of 
this classification and characterization of BF theory is, 
thus, phenomenological. Nevertheless, we will find the 
correct BF theory for the Z2 spin liquids with the gap- 
less edge states. More specifically, we will find that the 
BF theory matches the the crystal momenta and topo- 
logical degeneracies, and the nature of the edge theories 
of Z2 spin liquids. This justifies the correctness of the 
effective BF theories. 



(A) 



(C) 

E{k x ) 



(B) 



(D) 




E(k x ) 

> % 



FIG. 2. Stability of the gapless edge states. (A) To gap 
out the dangling Majorana fermions, the Majorana fermions 
should pair up and dimerize. However, it is clear from the 
figure that the dimerization doubles the unit cell along the 
boundary. This is similar to th e physics of surface states of 
weak topological insulator d 52 * 53 *. (B) Coupling of the dangling 
Majorana fermions to the quantum Ising chain at criticality. 
As the Hilbert space per site of the Ising chain has the di- 
mension 2 = (\/2) 2 , the Ising chain can be thought as the 
chains of two Majorana fermions per site. So when the Ising 
chain is coupled to the dangling Majorana fermions, there are 
three Majorana fermions per site which is guaranteed to be 
stable by counting Eq.(lll. (C) This picture can be further 
supported by looking at momentum space. The Ising chain 
at criticality is equivalent to a helical Majorana fermion state 
at k — tt and the edge state consist of the dangling Majorana 
fermions of the spin liquid has one mode at k = and k — tt. 
Thus there are total three modes at k = tt (k — is decoupled 
from k = 7r as the coupling between k — and k = n would 
double the unit cell). (D) When the interaction between the 
Ising chain and the edge state of the Z2 spin liquid is turned 
on, we will be left with one reconstructed Majorana mode at 
k — tt and another mode at k — 0. 



A. Crystal momenta and degeneracies of the Z2 
spin liquids 

To find BF theory for the Z2 spin liquids, we first need 
to compute the crystal momenta and degeneracies for the 
spin liquids. We assume that the mean-field Hamiltonian 
is enough for computing these quantities, i.e., the gauge 
fluctuation and interactions beyond the mean-field states 
do not change the crystal momenta and degeneracies. 

Here, we briefly show how to compute^ the crystal 
momenta and degeneracies of the Z2 spin liquids. The 
calculation is based on the mean-field theory for the Z2 
spin liquids. In the mean-field theory, we can assume 
that there are two fermions \I> T = (ip%,ip2) per site. To 
connect to the spin-1/2 Hamiltonian, we identify ipi = /j- 
and "02 = ft where t / 4- denotes the spin up and down 

in basis of S z , i.e., S = ^E a ,(34,j/^^//)- Tnis P ar " 
ton construction enlarges the Hilbert spaces {|^ me on >}> 



and we need to project it down to the physical states 
{I'i'spm >}• The projection is given by 

\®spin >= P\^rnean >, (15) 

where P projects the mean field state of fermions into 
the physical spin state with one fermion per site. In this 
paper, we specifically concentrate on the case of Zi spin 
liquid, where the fluctuating gauge field is Z 2 gauge field. 
According to projective symmetry group classificatiorPS 
there are two types of Z2 spin liquids if we have only 
translational symmetry: Z2A spin liquids (so-called zero 
flux states whose spinon band structure explicitly pre- 
serves translational symmetry) and Z2B spin liquids (so- 
called 7r-flux states whose spinon band structure doubles 
the unit cell). Here we focus on Z2A spin liquids. We 
introduce the fermion operator for 'S'k in the momentum 
space 'i'k = (^i,k, ip\ _k' Vfyki V4 -k) T ano - write down 
the Hamiltonian for ^k- Then, we see that ^J^k = 2 
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(number constraint) and ^l^cr ® <J 3 ^k = (particle- 
hole symmetry) for all k. In the presence of translational 
symmetry, phase transitions between different Z 2 spin 
liquids happen when there is a level crossing at the high- 
symmetry points of the Brillouin zone^l. There are four 
such points: (0, 0), (0, 7r), (tt, 0), (tt, tt). At these points, 
fermion occupation number = <S) o^^k can take 

different values at these time-reversal-invariant momenta 
(TRIM) while — 2 for all other possible k in the Bril- 
louin zone^. Different fillings at TRIM are the origin 
of the crystal momenta of Z2 spin liquids. The ground 
states can be obtained by filling all the states below the 
chemical potential fi = (half-filled). Hence the ground 
state carries the crystal momentum depending on the fill- 
ing rtk at these high-symmetry points. 

Naively, there are four ground states for a Z2 
spin liquid on the torus, and the four ground states 
can be labelled by the boundary conditions along 
the two directions of the torus (the four ground 
states are generated by imposing 'periodic' x 'periodic', 
'periodic' x 'anti-periodic', 'anti-periodic' x 'periodic', 
'anti-periodic' x 'anti-periodic' boundary conditions 
along the two directions of the torus to the fermions). 
As the total number of the fermions in the spin liquid 
should be even and thus we require X^kesz nk ^ 
to be a physical state. It is not difficult to see 
that the requirement X^kesz nk G 2Z should re- 
duce to Y^keTRLM n k S 2Z because = 2 for all 
ke BZ\ TRIM 19 . However, not all of k e TRIM is al- 
lowed for the fermion. For example, on the odd-by-even 
lattice (with the size of the system N x x N y such that 
N x € 2Z + I and N y € 2Z) with the periodic boundary 
conditions along the two directions of the torus, the 
allowed momenta (k x ,k y ) for the fermions can take the 
values of (^, ^) with n x = 0, 1, 2, • • • N x - 1, N x 
and n y = 0, 1, 2, • • • N y — l,N y . In this case, there is 
no allowed state for the fermions at k = (ir,0) and 
(tt, tt), and thus the summation X)keT_R/M n k should 
be modified to the summation X)keTi?/Af* nk wnere 
TRIM* is the set of the allowed momenta among 
k 6 TRIM. On the other hand, if the boundary con- 
dition along i-direction is changed to be anti-periodic 
while keeping the boundary condition along y - direction 
periodic, then the allowed momenta (k x ,k y ) is changed 



to k x 



2tv 
N, 



(n x + \),n x 



0,1,2,- --Nx with k y as 



before. For this case, TRIM* is identical to TRIM i.e., 
riirM* = {(0,0),(7r,0) ) (0,7r),(7r,7r)}. 

In the given lattice with the specificed boundary con- 
dition, we can compute the sum SkeTfl/M* nk ^ or a ^ 2 
spin liquid. If J2keTRiM' "-k £ 2Z + 1, then the state la- 
belled by the boundary condition is not a physical state 
and the state should not be considered as one of the 
ground states. Thus, the number of the ground states on 
the lattice can be smaller than the naive expectation 4 on 
torus. When YlkeTRiM* nk e ^Z, * nen ^ ne crvs t a l mo_ 
menta for the state is given by SkeT_R/M* krtk- Hence, 
if a Z2 spin liquid is given, we can compute the number 



of ground states on the four lattices (even-by-even, even- 
by-odd, odd-by-even, and odd-by-odd lattices) and label 
the ground states of the four lattice with the crystal mo- 
menta. For example, there are only two ground states for 
Wen's plaqucttc model on the even-by-odd lattices, and 
the two states carry the crystal momenta (0, 0) obtained 
by computing ^2keTRiM* ^ nk ^ or t ne ^ w0 ground states. 
We can write this compactly as follows (we denote "exo" 
lattice as even-by-odd lattice) 

(ex) Wen's plaquette model 
exo lattice: (0,0), (0,0) 

Following the above discussion, we find the crystal 
momenta of the ground states for the three Z2 spin 
liquids in the previous section (here we again denote 
"exe lattce" as even- by-even lattice, and "exo lattice" 
as even- by-odd lattices, etc.) 

(1) Wen's plaquette model, and n = x ±y 



exe lattice: 
exo lattice: 
oxc lattice: 
oxo lattice: 

(2) ft = x 

exe lattice: 
exo lattice: 
oxe lattice: 
oxo lattice: 

(3) n = y 

exe lattice: 
exo lattice: 
oxc lattice: 
oxo lattice: 



TT, 70,(0,0), (0,0), (0,0) 

0,0), (0,0) 
0,0), (0,0) 
0,0), (0,0) 

0, tt), (0,0), (0,0), (0,0) 
0,0), (0,0) 

0,^,(0,0), (0,0), (0,0) 
0,0), (0,0) 

tt,0), (0,0), (0,0), (0,0) 
^,0), (0,0), (0,0), (0,0) 
0,0), (0,0) 
0,0), (0,0) 



We will see that this complicated pattern of the crystal 
momentum can be reproduced by the effective BF theory 
in the next subsection. While deriving the above result, 
we implicitly assumed that the Hamiltonian Q({f}) = 
±1 for the trivial chain not to have any potential struc- 
ture in it. While the Hamiltonian Q here is artificial, we 
need to close the gap at the high-symmetry points of the 
Brillouin zone to change the crystal momenta. Hence, 
the pattern of the crystal momenta should be properties 
of a phase, not of a point in the phase diagram. 



B. Classification of BF theory on the square lattice 
with translational symmetry 

To classify BF t heory, we begin with the lesson from 
the previous studies^^ that we will utilize in this sec- 
tion. The first lesson is that encoding the transformation 
on the gauge theory is not 'gauge' degree of freedom, but 
it is a physical symmetry. For example in the toric code 
or Wen's plaquette model, we exchange the electric ex- 
citation and the magnetic excitation under the transla- 
tion. This is different from SX(2;Z) gauge symmetry for 



BF theory. The second lesson is that we can obtain the 
complicated patterns for the crystal momenta and degen- 
eracies by implementing the non-trivial transformations. 
Hence, the strategy is simple: finding all the possible 
physical transformations for BF field theory. In the fol- 
lowing discussion, we will systematically search and im- 
plement the transformations for BF theory. 

We begin with the Coulomb gauge ao = bo = for 
the gauge fields a M and 6 M in Eq. 0. Then, this reduces 
Eq.fJlJ into the following form (without sources) 



not physical) solution Q is 



Q 



which acts as (a x ,b x ) 



1 

1 



-> 







-1 






-1 





(22) 



(a y ,—b y ) and 



{a y ,b y ) 



L 



2tt 



v T ■ J ■ <9qv 



(16) 



(a x , —b x ). The fourth non-physical solution P is 

10 
0-100 
0-10 
1 



(23) 



where v T = (a x ,a y , b x ,b y ) and 



+: 
0-10 
0+100 
-10 



(17) 



which flips the sign of a y and b x spontaneously. A 
physical argument for excluding P and Q comes from 
the coupling of BF fields to the Z 2 charge field z and 
the Z2 vortex field v. We have the full Lagrangian 
L = Lbf + L rmm where Lbf is Eq.(16) and (in the 
= b = 0) 



"BF T J-ieoup 

Coulomb gague ao 



To be invariant under a unit lattice translation t x in x- 
axis (and similar for t y in y-axis), we require 



= \(d j -ia j )z\ 2 + \{d j -ib j )v\ i 



(24) 



L 



1 ^ 

27T 1 



•j-a v = i-(t, [v T ]).j-a (^ [v T ]). (is) 



Now, we translate the system by one unit lattice along 
x-axis which acts as t x on the fields. Here, t x [d^] = d t ,. 
and t x [Lbf] — Lbf , and we are required to satisfy the 



Hence, we look for solutions t x and t y satisfying Eq.([l8]). equation L. 



trr- \L r 



where 



The obvious possibility for t x and t y is to consider the lin- 
ear transformations acting on v, i.e., we associate (t x ,t y ) 
with the matrices (O x , O y ) and the vectors (u^, u y ) such 
as 



\L C 



\{dj - it x [aj])t x [z] | 2 + \(dj - it x [bj])t x [v] 



(25) 



: v — > v. For t x = 
z. However, for P 



(19) 



coup\ 



We plug this into Eq.(18) and find a restriction on O x 



J. 



(20) 



The restriction for u x , 9 shows up when we take the mi- so should be u^. More specifically, [ 



croscopic picture into consideration. Before proceed- 
ing further, we notice that XiV is orthogonal, i.e., 
det(O a> j / )=det(Oj ) = ±1 as J is invertible. This alone 
cannot fix the form of O x y . We also note that O x — 1. 
This is because the size of the smallest even-by-even lat- 
tice is 2x2, and we expect the degeneracies on the even- 
by-even lattice are trivially 4 i.e., we force the gauge field 
v to return to itself. Then, it is easy to check that there 



For t x = I, we have t x : z — > z and t a 
O*, we have t x : z — > v and t x : v — > 
and Q, there is no way to make the equality t x [L, 
L C oup from t x [z] and t x [v] , hence we exclude them as the 
physical transformation of BF theory which is consistent 
with the Z2 spin liquids in the previous section. 

Now, we move on to the constant vector u XyV . As the 
vectors v are Z 2 gauge theory in the microscopic picture, 

{0,^} 



u *,yla € 



are four solutions for O x _ y satisfying Eq.(20) (Note that 



the overall sign of the O x ^ y is irrelevant as the sign of Z 2 
theory is not important). The first physical solution is 
the trivial I =diag(l, 1, 1, 1). The second physical solu- 
tion is 'twist' O* 



mod 2tt. However, not all of u x ^ y pattern are physically 
independent. There are n flux or no flux in Z 2 gauge 
theory for and bj, and this reduces many u XtV patterns 
into only four of them. Those are the followings: 1) 
Both and b v have fluxes in the unit cell. 2) only a M 
contains fluxes in the unit cell. 3) only b v contains fluxes 
in the unit cell. 4) No fluxes are present. For 1), we have 
" = (0,0,tt,0). For 2), we have 
\ = (0,0,0,0). For 3), we have 
(0,0,tt,0). For 4), we have 



ui = (0, 7T, 0, 0) and u 



u x = (0, 7r, 0, 0) and u. 
T = (0,0,0,0) and 



u 



u 



x 



V 



u: = (0,0,0,0). 



o' 



10 
1 
10 
10 



(21) 



—¥ bj and 6„ — > 



which exchanges and &j, i.e. 
under a unit lattice translation. The third possible (but 



We write (t x ,t y ) as {(O x , u x ), (O y , u y )}. We denote 
twist matrix simply as O* and trivial matrix as I, and 
we also denote the nontrivial flux vector as u and no flux 
as 0. For example, we can have BF theory with the 
translational invaraint property {(O*, 0), (I, Uy)} satis- 
fies the following relations: t x : (cii,bj) — > {bi,aj) and 
t y : (a i} bj) -> (ai,bj +tt5 x j). 

A simple combinatorial computation then concludes 
that this classification gives 2x2x4= 16 classes for BF 
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theory consistent with the underlying Z2 gauge theory 
on the square. Among 16 classes, six classes are not of 
interest. To see this, we look at the commutation relation 
for t x and t y 

[t x ,t y }^0 (26) 

seemingly implying that those states are not physical. 
But this is not true; for example, the magnetic transla- 
tions are not commuting each other in the quantum Hall 
states. Still, there is no allowed ground stated in the 
odd- by-odd lattices if [t x ,t y ] ^ 0. As all of our targeting 
states have 2 or 4 states, we exclude the BF theory of 

[t X ,ty] ^0. 



C. Crystal momenta, degeneracies, and Z2 indices 

Upon obtaining the transformations for the relevant 
BF theory, we compute the crystal momenta of the 
ground states on even-by-even, odd-by-even, even-by- 
odd, and odd-by-odd lattices. We follow the straightfor- 
ward calculation in the reference [T5] to get the crystal 
momenta. 

Before presenting a series of results from the compu- 
tation, we exhibit two examples showing how we get the 
crystal momenta on lattices. 

Example 1. {(I, u x ), (1, 0)} 

We have the transformations under the translations as 
t x ■ (ai,bj) -> (oj +n& y ,i,bj) and t y : (04,67) -> (cn,bj). 
We denote the zero mode of BF fields (ai,bj) as 
(6i,4>j). Then, the BF theory relates 9 X (<j) x ) to <p y 
{By) as the canonical conjugate pairs. Explicitly, we 
have [0 x ,4>y\ = iit and [4> x ,0 y ] = i7r. Then, we have 
four well-defined ground states (|1 >,|2 >,|3 >,|4 >) 
such as |2 >= e- i9 "\l >, |3 >= e _ ^*|l >, |4 >= 
e -te„ e -%4>^i > wn ere |1 > transforms trivially un- 
der any transformations of t x , y - Then, we have 
t x : (|1 >, |2 >, |3 >, |4 >) -> (|3 >, |4 >, |1 >, |2 >) and 
t y : (|1 >,|2 >,|3 >,|4 >) (|1 >,|2 >,|3 >, |4 >). 
Hence, [t x ,t y ] =0 and t x and t y can be simultaneously 
diagonalized. We conclude that the two states of the 
form ^=(|1 > -|3 >) and ^=(|2 > -|4 >) carry the 

crystal momentum (tt, 0), and other two states carry the 
crystal momentum (0, 0) for any lattice. This gives the 
crystal momentum spectrums (n, 0), (n, 0), (0, 0), (0, 0) 
on even-by-even, even-by-odd, odd-by-even, and odd-by- 
odd lattices, and there is no matched Z2 spin liquids of 
interest. 

Example 2. {(O*,0), (1,0)} 

We have the transformation law under the translations 
t x : (di,bj) -> (6j,Oj) and t y : ((H,bj) -> ((H,bj). We 
begin with the even-by-even lattice. Then, we can 
expand BF theory and obtain the zero modes 6 X (4> x ), 
canonical conjugate to (f) y (9 y ). Then as before, we have 
four well-defined ground states (|1 >,|2 >,|3 >,|4 >) 



such as |2 >= e~ ie -\l >, |3 >= e -^|l >, |4 >= 
e -i0x e -i<l>x\\ > where |1 > transforms trivially un- 
der any transformations of t xy . Then, we have 
t x :(\l >, |2 >, |3 >, |4 >) -> (|1 >, |3 >, |2 >, |4 >) and 
t y : (|1 >,|2 >,|3 >,|4 >) -> (|1 >,|2 >,|3 >, |4 >). 
Hence, [t Xl t y ] —0 and t a and t y can be simultaneously 
diagonalized. We conclude that one state of the form 
^(|2 > — 13 >) carries the crystal momentum (tt, 0), 
and the other three states carry the crystal momentum 
(0, 0) for even-by-even lattice. This gives the spectrum 
of the crystal momenta (71", 0), (0, 0), (0, 0), (0, 0) on 
even-by-even lattice. Now, we move to the even-by-odd 
lattice. In this case, we can expand the zero mode 
as well as in even-by-even lattice because t y is triv- 
ial, and this gives the crystal momentum spectrum 
(tt,0), (0,0), (0,0), (0,0) on even-by-odd lattice. It 
turns out to be different for odd-by-even lattice as the 
non-trivial boundary condition T x : (eii,6,) — > (6j,£jj) 
where T x = (t x ) LlE , and we cannot expand BF theory 
to obtain the zero mode. The resolution for this is to 
double the odd-by-even lattice along x- axis to get even- 
by-even lattice. Then, we can now expand BF fields to 
obtain the zero modes (|1 >, |2 >, |3 >, |4 >) such as 
|2 >= e~ ie *\l >, |3 >= e- 4 ^|l >, |4 >= e - ie -e- l ^\l > 
with the boundary condition T x : (0 x ,(f> x ) — > (4> X ,0 X ). 
The boundary condition implies that the originally 
four independent states bind to each other to form 
only two independet states |1 >, |2 >. Furthermore, 
t x : (|1 >,|2 >) (|1 >, |2 >), and thus the crystal 
momentum spectrum on the odd-by-even lattice is 
(0,0), (0,0). The similar consideration on odd-by- 
odd lattice gives the crystal momentum spectrum as 
(0,0), (0,0). These crystal momentum spectrums match 
the spin liquids of h — y in the previous section. 



This calculation can be also done for other transla- 
tional symmetric BF theory, and now we find that 

(1) ({(O*,0),(O*,0)}) is Wen's model or n = x ± y 

(2) ({(O',0), (1,0)}) is h = yZ2 spin liquid 
(2) ({(I, 0), (O*, 0)}) is n = x Z2 spin liquid 

It is not difficult to see that these BF theories have 
the same edge spectrum as the micrscopic consideration 
in the previous section. The spectrum of the crystal mo- 
menta and degeneracies for other BF theories can be 
found in the appendix. 

Before finishing this section, we comment on the con- 
nection to the previous studjEil. The pattern of the crys- 
tal momenta and degeneracies of Z2 spin liquids are, in 
fact, fixed by the fermion parity at the high symmetry 
points in the Brillouin zone, and this fact can be used to 
classify Z2 spin liquids. We can find that Wen's model 
and ii = x ±y belong to [0110] and [1001] classes in the 
classification of Z2 spin liquids in the index system of 
the previous study (h — x belongs to [1100] and [0011] 
classes, and h = y belongs to [1010] and [0101] classes.) 
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Though it can be shown that these Z2 spin liquids with 
the gapless edge belong to the particular class, it is not 
clear if this index system implies the gapless edge states 
between Z2 spin liquids and vacuum. It would be an 
interesting future direction to study if the indices imply 
the presence of the gapless edge state. 



V. GENERALIZATIONS AND CONCLUSIONS 

Can we generalize our reasoning for the gapless edge 
states of Z2 spin liquids on the square lattice to other lat- 
tices such as the triangular lattices? More interestingly, 
can we generalize it to higher dimensions? It is not diffi- 
cult to note that the answers to these questions are 'y es '> 
at least in the mean-field Hamiltonian. The nature of the 
edge states will depend on the direction of stacking one- 
dimensional Majorana fermion chains for f,^ u and fi^d- 
Hence, we can apply this to any two-dimensional lattice 
to construct Z2 spin liquids with the gapless edge states. 
We can also generalize to three-dimensional spin liquids 
to write down the mean field Hamiltonian, but it is bet- 
ter to have an exactly solvable model that exhibits the 
gapless surface states for the cubic lattice. The model is 
the direct generalization of Wen's plaquette model on the 
square lattice, and it involves six Majorana fermions per 
site, i.e., this is a Hamiltonian for a spin-3/2 system 5 -^. 
We start with the I2 gauge theory on the cubic lattice 



H 



-<] 



e^=-*e n u * 



(27) 



P <ij>GP 



where P denotes the plaquette of the cubic lattice. Here, 
Uij = ±1 is the Z2 gauge on the link < ij >, and 
thus n<y >ep Uij — Fp is the field strength on the pla- 
quette P. The particular (but general) representation 



of this gauge theory is U, 



i.i+x 



i\ ity Xi+y,y, Ui^+z = i\i, z Xi+z,z with the translational 
symmetry. Hence, we have the six Majorana fermions 
(\,x, K, v , K,z,Xi,x,Xi,y,Xi,z) per site, and the dimension 
of the Hilbert space per site is 2 6 / 2 = 8 states per site. 
Again we need to halve this Hilbert space to represent a 
spin-3/2 system,, i.e., the Hilbert space per site should 
be 4. This is done by taking the Z2 redundancy into ac- 



count: Eq.(27) is invariant under Uij — > SiUijSj, Sj = ±1 



In the three-dimensional model, we see that all the Ma- 
jorana fermions in the bulk are paired. On the boundary, 
however we have one dangling Majorana fermion per unit 
cell on the boundary. Following the reasoning in the two- 
dimensional system, the flat bands on the boundary will 
start to disperse to form the gapless surface states when a 
small perturbation away from this exactly solvable limit 
induces 



P^p.dae. — 



(txViVi+x + tyViVi+v)- 



(28) 



This Hamiltonian has the gapless spectrum E(]s.) = 
Jt 2 sm 2 (k x ) + t 2 sm 2 (ky). To gap out the spectrum, we 



need to double the unit cel l 52 * 53 * which is prohibited by 
the translational symmetry. 

As before, we can deform the exactly solvable Hamilto- 
nian to obtain a series of anisotropic Z2 spin liquids with 
the gapless surface states. How many different classes 
of the Z2A spin liquids with translational symmetry in 
the cubic lattice are there (if two Z2 spin liquids have 
the same crystal momentums and degeneracies, and the 
surface states, we define them as the same class)? We re- 
peat the argument in the two-dimensional case and con- 
struct 3D mean-field (ideal) Hamiltonians of Z 2 spin liq- 
uids from the one-dimensional Majorana fermion chains. 
By checking the crystal momenta and the degeneracies 
in the three-dimensional lattices, we found that there are 
14 classes which have the gapless surface states protected 
by the translational symmetry. The effective theory for 
these Z2 spin liquids should be (3+1)- dimensional BF 
theor^ 



L 



1 

27 



2^ °v 



(29) 



which is known to have the gapless Dirac-like surface 
spectrum^. It would be an interesting future research 
direction to show how the translational invariance stabi- 
lizes the gapless surface states of (3+l)-dimcnsional BF 
theory . 

In summary, we considered translationally-symmctric 
Z2 spin liquids which have gapless edge/surface states. 
The edge/surface states are constructed out of the dan- 
gling Majorana fermions on the boundary. We clarified 
the conditions for the gapless edge states for Z2 spin liq- 
uids and constructed the effective BF theory reflecting 
the underlying translational symmetry of the lattice. The 
general reasoning applies to the three-dimensional lat- 
tices, and we demonstrated that there is a gapless surface 
state. 
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CRYSTAL MOMENTA OF BF THEORY 

We list the crystal momenta of each theory (we denote 
"exe lattce" as even- by-even lattice, and "exo lattice" 
as even-by-odd lattices, etc) 

(1) ({(1,0), (1,0)}) is [0000] Z2 spin liquid 

(2) ({(I,u x ),(I,0)}) 
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exe lattice: (tt, 0), (tt, 0), (0, 0), (0, 0) 
exo lattice: (tt, 0), (tt, 0), (0, 0), (0, 0) 
oxe lattice: (tt, 0), (tt, 0), (0, 0), (0, 0) 
oxo lattice: (tt, 0), (tt, 0), (0, 0), (0, 0) 

(3) ({(I,0),(I,u y )}) 

exe lattice: (0, tt), (0, tt), (0, 0), (0, 0) 
exo lattice: (0, tt), (0, tt), (0, 0), (0, 0) 
oxe lattice: (0, tt), (0, tt), (0, 0), (0, 0) 
oxo lattice: (0, tt), (0, tt), (0, 0), (0, 0) 

(4) ({(I,u x ),(I,u y )}) is [1111] Z2 spin liquid 

(5) ({(O*, 0), (O*, O)}) is [0110] , [1001] Z2 spin liquid 

(6) ({(O*,iO,(O',0)}) ([t x ,t v ]?0) 

(7) ({(OSO),^*,!!,)}) ([*„*„] ^0) 

(8) ({(0*,u x ),(0*,u y )}) 

exe lattice: (tt, it), (tt, 0), (0, tt), (0, 0) 
exo lattice: (0,0), (0,0) 
oxe lattice: (0,0), (0,0) 
oxo lattice: (0,0), (0,0) 

(9) ({(O*,0), (1,0)}) is [0101] , [1010] Z2 spin liquid 



(10) ({(()', O), (I, u„)}) ([t x ,t y ]^0) 

(11) ({(()*, u*), (1,0)}) 

exe lattice: (tt, 0), (tt, 0), (0, 0), (0, 0) 
exo lattice: (tt, 0), (tt, 0), (0, 0)(0, 0) 
oxe lattice: (0,0), (0,0) 
oxo lattice: (0,0), (0,0) 

(12) ({(0',iO,(I,u y )}) ([t x ,t y }^0) 

(13) ({(I,0),(O*,0)}) is [1100] , [0011] Z2 spin liquid 

(14) ({(1,10, (O*,0)}) ([«»,<»] ^0) 

(15) ({(I,0),(O*, U!/ )}) 

exe lattice: (0, tt), (0, tt), (0, 0), (0, 0) 
exo lattice: (0,0)(0,0) 
oxe lattice: (0, tt), (0, tt), (0, 0), (0, 0) 
oxo lattice: (0,0), (0,0) 

(16) ({(I.u^^O*,^)}) ([t x ,t y }^0) 

We looked at the patterns of the crystal momenta of 
the Z2 spin liquids^ and found four Z2 indices of Z2 spin 
liquids which match these patterns. Note that the classes 
(5) — (16) of BF theories have the gapless edge state. 
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